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Nomenclature 


Abstract 

Flowfield solutions over the Mars Pathfinder 
Probe spanning the trajectory through the Martian 
atmosphere at angles of attack from 0 to 11 degrees 
are obtained. Aerodynamic coefficients derived from 
these solutions reveal two regions where the derivative 
of pitching moment with respect to angle of attack is 
positive at small angles of attack. The behavior is asso- 
ciated with the transition of the sonic line location be- 
tween the blunted nose and the windside shoulder of the 
70 degree half-angle cone in a gas with a low effective 
ratio of specific heats. The transition first occurs as the 
shock layer gas chemistry evolves from highly nonequi- 
librium to near equilibrium, above approximately 6.5 
km/s and 40 km altitude, causing the effective specific 
heat ratio to decrease. The transition next occurs in an 
equilibrium flow regime as velocities decrease through 
3.5 km/s and the specific heat ratio increases again with 
decreasing enthalpy. The effects of the expansion over 
the shoulder into the wake are more strongly felt on 
the fustrum when the sonic line sits on the shoulder. 
The transition also produces a counter-intuitive trend 
in which windside heating levels decrease with increas- 
ing angle of attack resulting from an increase in the ef- 
fective radius of curvature. Six-degree-of-freedom tra- 
jectory analyses utilizing the computed aerodynamic 
coefficients predict a moderate, 3 to 4 degree increase 
in total angle of attack as the probe, spinning at ap- 
proximately 2 revolutions per minute, passes through 
these regions. 
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C'a axial force coef. 
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D probe diameter, m 

h altitude, km 

M Mach number 

M mixture molec. weight, kg/kg-mole 

p pressure, N/m 2 

q convective heating, W/cm 2 

Q dynamic pressure, N/m 2 

r radius from probe axis, m 

Re Reynolds number 
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a angle of attack, degrees 
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Introduction 


The Mars Environmental Survey (MEStIR) 1 
Pathfinder 2 program proposes to observe the Martian 
atmosphere and surface with unmanned probes. Re- 
cently renamed Mars Pathfinder, the program, embrac- 
ing the design principles “faster, better, cheaper”, now 
focuses on entry, descent and landing of a single probe. 
The design program naturally borrows much, includ- 
ing body shape, from the earlier Viking project 3 (and 
associated ground-based experimental data 4- '), which 
landed a science station on the Martian surface in 1978. 
Unlike Viking, Pathfinder seeks to simplify the mission 
by utilizing direct, entry into the Martian atmosphere 
and eliminating any reactive control system during its 
descent. The direct entry forces a higher entry velocity 
and heating rates than experienced by Viking. Aero- 
dynamic stability is to be maintained at a nominal zero 
degree angle of attack with an on-axis cent.er-of-gra.vit.y 
location while the probe spins at approximately two 
revolutions per minute. 

While much of the research in support of the 
Viking project can be used for Pathfinder design val- 
idation, the differences in entry velocity and angle 
of attack necessitate a re-evaluation of Viking shape 
aerothermodynamics at Pathfinder conditions. Viking 
had been designed to fly at an eleven degree angle of 
attack, a condition for which the aerodynamic uncer- 
tainties associated with gas models at hypersonic con- 
ditions were relatively small for that shape. The pri- 
mary goal of this paper is to provide understanding of 
the aerodynamics for the Viking shape near zero de- 
grees angle of attack for the Pathfinder mission. It 
also demonstrates how heating distributions change as 
a function of angle of attack, and highlights the behav- 
ior of both the aerodynamic coefficients and heating 
distributions as a function of the sonic line location 
on the probe. The heating analysis contained herein 
is purposefully kept very simple; more comprehensive 
heating analyses which ignore the three-dimensional is- 
sues treated here may be found in the literature. 8-15 

To this end, a matrix of fast running, inviscid, per- 
fect gas solutions is generated and aerodynamic coef- 
ficients are extracted. The matrix spans angles of at- 
tack from zero to eleven degrees and trajectory points 
from peak dynamic pressure to parachute deployment. 
The matrix includes velocities for which ballistic range 
tests indicated static stability problems with the Viking 
shape at small angles of attack. The perfect gas, con- 
stant. specific heat ratio is chosen to match the equi- 
librium, normal shock density ratio at the specific 
freestream velocities and temperatures. Viscous, real 
gas solutions are concurrently generated on a comple- 
mentary matrix to: ( 1 ) ascertain the sensitivity of blunt 
body aerodynamics to the simplified perfect gas solu- 
tions; (2) define probe aerodynamics in the nonequi- 



Figure 1: Diagram of the Mars Pathfinder probe. 

librium flow regime prior to peak heating; and (3) es- 
tablish the convective heating distributions. Limited 
validation of results is provided with grid convergence 
studies and comparison to experimental data. 


Configuration and Trajectory 

The Mars Pathfinder configuration is shown in Fig. 
1. The forebody is a 70° half angle cone with base di- 
ameter 2.65 m, nose radius 0.6625m and shoulder ra- 
dius 0.06625 m. The on-axis center-of-gravity is located 
0.662 m behind the nose. Only the forebody is included 
in this study unless otherwise specified. ThgTorebody 
analysis is sufficient, for computation of aerodynamic 
coefficients at. hypersonic conditions. Ba.seflow effects 
are important, to computation of aerodynamic coeffi- 
cients at. low supersonic conditions and will be briefly 
discussed. 

A preliminary t.ra.ject.ory lb is generated assuming a. 
mass of 418 kg. and a. 7.6 krn/s ballistic entry at - 
18 degrees. Test, points along this trajectory, defined 
in Table 1, range from 18 s prior to peak heating to 
parachute deployment. 


Algorithms 

Perfect. Gas, Inviscid Analysis 

The perfect, gas, inviscid analysis is provided by 
the High ALpha. Inviscid Solution (HALIS) 1 ' code. 
The code employs an explicit. Ma.cCorma.ck predictor- 
corrector scheme. The bow shock defines the inflow 
boundary, on which the Ra.nkine Hugoniot. relations for 
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a moving shock define the inflow conditions. The invis- 
cid wall boundary condition enforces the tangency con- 
dition on the surface and corrects velocities to enforce a 
constant total enthalpy. Backward differences are used 
on the outflow boundary, which lies downstream of the 
sonic line in a supersonic domain. High confidence in 
this code for inviscid, blunt body flow computations is 
derived from extensive comparisons as documented in 
the literature. 18-20 

The perfect gas analysis requires an approxima- 
tion to real-gas chemistry behind the bow shock using 
an appropriate choice for the ratio of specific heats, 7. 
A standard approximation for air chemistry based on a 
local ratio of enthalpy to internal energy 21 causes some 
problems when used for Martian atmospheric gases 
because of the large, negative heat of formation for 
CO2. An alternative approach used here is to choose 
a constant value for 7 which matches the density ratio 
(rj = P 2 I Pi) across a normal shock for equilibrium 
flow. 

7 = 2 -H ( 1 ) 

V ~ 1 (r, - 1 )Ml V ; 


The values for rj were computed using a viscous shock 
layer code 12 with free energy minimization for the Mar- 
tian atmospheric chemistry. This approximation is only 
applied to points on the trajectory at or following peak 
dynamic pressure where the equilibrium normal shock 
density ratio was in good agreement with the nonequi- 
librium calculation discussed below. These effective 
values for 7 in the equilibrium regime are presented 
in Table 1. Sound speed is calculated with the perfect 
gas relation: 


Real Gas, Viscous Analysis 

The real gas, viscous analysis is provided by the 
Langley Aerothermodynamic Upwind Relaxation Al- 
gorithm (LAURA). 22-24 Comparisons to experimen- 
tal data for hypersonic flows in air are documented in 
the literature. 25-29 The code employs upwind-biased, 
point-implicit relaxation. Inviscid fluxes are approxi- 
mated with Roe’s averaging, 30 eigenvalue limiting (sim- 
ilar to Harten 31 ), and Yee’s symmetric total variation 
diminishing scheme. 32 Viscous fluxes are approximated 
with central differences. A recently added model for 
Martian atmospheric constituents includes 8 species 
(CO2, CO, N 2 , O2, NO, C, N, O). 14 Reacting gas 
chemistry is used for the hypersonic cases (Mach 9 and 
above). The code is sufficiently robust to achieve con- 
vergence at near equilibrium conditions. Frozen chem- 
istry is employed for the Mach 1.9 case because the post 
shock conditions are not severe enough to induce chem- 
ical reactions. Thermal equilibrium is specified for all 
points below peak dynamic pressure (Mach 22.3). A 
two-temperature model is applied at trajectory points 


prior to and including peak dynamic pressure. The 
local Mach number is based on a frozen (species, vibra- 
tion, electronic excitation) sound speed: 33 


a 
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Unlike HALIS, LAURA employs shock capturing 
and the inflow boundary and grid are aligned with the 
captured bow shock. No-slip conditions are applied 
in a pseudo cell behind the surface within the context 
of a finite-volume formulation. Wall temperatures are 
held constant and defined in Table 1. Species mass 
fractions at the wall are held fixed at their respective 
freestream values ( cco 2 = -97, cjy 2 = -03, all others 
set to 1CT 10 ). Extrapolation is used to close the outflow 
boundary. 

LAURA can be used to generate inviscid, perfect 
gas solutions. However, the wall boundary conditions 
in this mode are very simple, specifying zero normal 
mass flux and extrapolating pressure and energy to the 
wall with no corrections for total enthalpy conservation. 
Nevertheless, LAURA is used in an inviscid, perfect gas 
mode to ascertain sensitivity of the solution for aero- 
dynamic coefficients to algorithm at selected trajectory 
points. The differences are small and will be discussed. 


Computational Grid 

The HALIS solutions are generated on a cylin- 
drically derived coordinate system with the axis co- 
incident with the vehicle axis. The surface grid (Fig. 
2a) uses 46 points in the radial (predominantly stream- 
wise) direction and 36 points in the circumferential di- 
rection. The volume grid is built up in a body normal 
direction using 21 points between the surface and the 
bow shock. 

The surface-grid distribution for the LAURA test 
cases (Fig. 2b) derives from a singularity free mapping 
that provides a nearly uniform mesh in the stagnation 
region and enhanced circumferential and axial resolu- 
tion at the shoulders of generic aerobrake shapes. The 
surface grid may be perceived as a rectangular sheet 
which has one boundary on the symmetry plane and 
the other three boundaries stretched along the outer 
edge of the vehicle. Grid lines in the vicinity of the 
two sheet corners on the circular shoulder are some- 
what skewed. The baseline grid is defined by 61 x 31 
point rectangle which yields 61 points along the symme- 
try plane and 121 points along the outer, semi-circular 
edge. The volume grid is built up in a body normal 
direction using 65 points between the surface and the 
inflow boundary which lies upstream of the captured 
bow shock. The normal distribution is defined with a 
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Figure 2: Forebody surface grid on the Mars Pathfinder 
probe. 

cell Reynolds number equal to 1 at the surface and a 
stretching function in which roughly half the available 
cells are assigned to the boundary layer and the remain- 
der are equally spaced across the shock layer. 34 Normal 
grid distribution is periodically updated as the solution 
evolves. Some grid refinement studies are implemented 
using an 81 x 41 point surface grid. 

Differences in thgttreatment of the bow shock pre- 
clude use of the same grid for both solution techniques. 
Both grid structures emphasize resolution in the vicin- 
ity of the circular shoulder because details of the ex- 
pansion process are critical to the evaluation of aero- 
dynamic coefficients. The rectangular topology, though 
more difficult to generate, is superior to the cylindri- 
cal topology for the purpose of resolving pressure and 
heating distributions around the vehicle axis, unencum- 
bered by the presence of a coordinate singularity. 


Aerodynamics 

HALIS: Inviscid, Perfect Gas Aerodynamics 

The inviscid, perfect-gas solution matrix run with 
the HALIS code and the associated aerodynamic coef- 
ficients are presented in Table 2. The matrix includes 
angles of attack of 0° (design condition), 2°, 5°, 8°, and 
11° (Viking condition). A base pressure equal to p m 
is assumed in the computation of aerodynamic coeffi- 
cients. 

Lift, drag, and pitching moment coefficients are 
plotted as a function of Mach number in Fig. 3 with 
the solid lines and filled symbols. Recall that the ef- 
fective specific heat ratio 7 also varies with Mach num- 
ber, so that gas chemistry effects are approximately 
included in these results. The lift coefficient is insensi- 
tive to Mach number for this configuration above Mach 


6. The drag coefficient shows a stronger dependence on 
Mach number. The pitching moment coefficient shows 
a distinct hump (local maximum) as a function of Mach 
number and the location of the hump changes with an- 
gle of attack. 

LAURA: Viscous, Reacting Aerodynamics 

The local maxima noted in the HALIS results sug- 
gests the possibility that the the pitching moment may 
be positive at small angles of attack near Mach 12 
in which case the vehicle would be statically unsta- 
ble. Concern regarding this feature prompted an ex- 
panded test matrix for the viscous, real-gas solutions 
using LAURA. The test conditions and associated aero- 
dynamics are presented in Table 3. Lift, drag, and 
pitching moment coefficients are plotted as a function of 
Mach number in Fig. 3 with the hollow symbols. These 
results include viscous forces. Aerodynamic coefficient 
magnitudes including unsteady baseflow effects (to be 
discussed subsequently) at Mach 1.9 and 5 degrees an- 
gle of attack are indicated by the open symbols within 
the bracket labelled base effects. The LAURA and 
HALIS predictions for lift coefficient are in excellent 
agreement for all Mach numbers and angles of attack 
tested. The LAURA and HALIS predictions for drag 
and pitching moment coefficient begin to diverge above 
Mach 9. The maximum percent difference in drag co- 
efficient, 1.3%), occurs in the test matrix at Mach 12, 
a = 2°. The maximum difference in pitching moment 
coefficient, 0.00245, occurs at Mach 16, a = 2°. The 
LAURA results predict static instability at this Mach 
number for small angles of attack. 

The ratio of pitching moment to normal coefficient 
is related to a nondimensional stability margin. This 
ratio, plotted in Fig. 4, collapses for all angles of at- 
tack in the test matrix below Mach 10 into a straight 
line. (Baseflow effects at Mach 1.9 are not included in 
these curves.) Above Mach 20, the angle of attack de- 
pendence diminishes again (until nonequilibrium effects 
become important as discussed in the next section). In 
between, mechanisms operate which are clearly sensi- 
tive to both Mach number (gas chemistry) and angle 
of attack. 

Pressure and Moment Distributions 

At small angles of attack, both C m and Cn are 
small numbers whose magnitude depend on the differ- 
ence of large windside and leeside forces. (The term 
leeside is somewhat of a misnomer; it refers to the side 
of the probe pitched away from the oncoming flow even 
though direct impingement still occurs.) To better un- 
derstand how these windside and leeside forces evolve 
to produce the static instability at a = 2°, pressure 
distributions on the symmetry plane are presented in 
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Figure 4: Ratio of pitching moment coefficient to nor- 
( a ) Drag coefficient, m al force coefficient, versus Mach number at, several 

angles of attack a for predominantly equilibrium flow 
regime. 
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Figure 3: Aerodynamic coefficients versus Mach num- 
ber at, several angles of attack a for predominantly equi- 
librium flow regime below M m = 22.3, maximum 

dynamic pressure. 


Figure 5: Pressure distribution along the symmetry 
plane for three Mach numbers spanning the static in- 
stability at, a = 2° in the predominantly equilibrium 
gas chemistry regime. 


Fig. 5 at trajectory points before the static instabil- 
ity (Mco = 22.3), at, the peak of the static insta- 
bility (Mco = 16.0), and after the static instability 

(Mco = 9.4). On the leeside (x > 0.5) the pressure 
for Mco = 22.3 has a, relatively constant,, positive slope 
all the way to the circular shoulder (x = 1.28) at which 
point, it, drops sharply in the expansion to the wake. 
The nondimensional pressure at, M m = 16.0 is approx- 
imately 1.5 % higher than at, M m = 22.3 but, shows an 
earlier and more moderate expansion. At, M m = 9.4 
the pressure distribution has a, stronger, negative slope 
and the pressure drop at, the shoulder feeds far up- 
stream. On the windside (x < —0.5) the pressure 

distribution starts relatively flat at, M m = 22.3 with 
some upstream influence of the shoulder evident, in the 
rounded drop at, x = —1.28. The flat distribution 

becomes more rounded, and average windside nondi- 
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mensional pressure levels continue to decrease with de- 
creasing Mach number as the effects of the expansion 
feed farther upstream. 

Though not easily seen, there is a small region near 
the shoulder where the leeside pressure exceeds the cor- 
responding windside pressure for the two highest Mach 
numbers. The behavior of the pressure distribution in 
this region significantly influences the pitching moment 
coefficient because of the relatively larger moment arms 
and surface area at the edge of the probe as compared 
to the inboard nose and fustrum regions. The inte- 
grated effect of pressure on C m at a = 2° is ap- 

proximately illustrated in Fig. 6 in which the product 
of pressure and moment arm (radius) is plotted as a 
function of area on the windside and leeside symme- 
try plane for M m = 22.3, 16.0, and 9.4, respectively. 
The integrated area under each curve, f prdA, repre- 
sents the pitch up moment (destabilizing) caused by 
leeside pressures and the pitch down moment (stabi- 
lizing) caused by windside pressures. The difference 
in the areas represents the net pitching moment. In 
all cases, the small value of a dictates that the differ- 
ence of the integrated areas will be small. However, 
as noted earlier, a crossover point occurs in the curves 
at the two highest Mach numbers (Figs. 6 a and 6 b). 
At Mco = 16.0 (Fig. 6 b) the crossover occurs far 

enough inboard such that tint leeside integrated area 
exceeds the windside integrated area, resulting in a net 
positive pitching moment. 

Influence of Sonic Line Location and 7 

Direct influence of the expansion over the shoulder 
on the cone fustrum pressure levels further upstream is 
regulated by the magnitude of the local Mach number, 
and limited completely to the subsonic zone bounded 
by the sonic line. Acoustic waves cannot travel up- 
stream in a supersonic flow. An exploded view of the 
sonic line location over the cone fustrum and circular 
shoulder is presented in Fig. 7a (leeside) and in Fig. 
7b (windside). At M m = 22.3 the leeside sonic line 
lies along the boundary layer edge, essentially elimi- 
nating upstream influence of the expansion, consistent 
with observations in Fig. 5. At M m = 16.0 the lee- 
side sonic line moves further out into the inviscid flow. 
A secondary subsonic bubble forms above the first, in 
the inviscid layer, its upper boundary embedded in the 
captured shock. At M m = 9.4 the leeside sonic line 
terminates in a line nearly perpendicular to the shoul- 
der, providing ample communication of the expansion 
to the entire shock layer as was evident in Fig. 5. 

The windside sonic lines in Fig. 7b all terminate 
at lines above the shoulder which curve sharply down- 
stream into the boundary layer in the rapidly expand- 
ing flow. The upper boundary of the windside sonic 
lines are embedded in the captured shock, and indicate 



(a) \l = 22.3 



(b) \l = 16.0 



(c) \l = 9.4 


Figure 6 : Approximation to the pitching moment coef- 
ficient integral in the vicinity of the leeside and wind- 
side symmetry planes at a = 2 °. 
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(a) Leeside. 



the significant influence of gas chemistry on shock layer 
thickness across this Mach number range. The size of 
the windside subsonic zone at M m = 9.4 is over twice 
as large as the subsonic zone at M m = 22.3, and 

the windside pressures in Fig. 5 again confirm that the 
larger subsonic zone yields larger pressure relief over 
the windside cone fustrum. 

The HALIS (perfect gas) predictions for sonic line 
location are not shown here, but they exhibit essen- 
tially the same behavior and relationship with pressure 
distribution as predicted by LAURA. Furthermore, the 
HALIS solutions are in reasonable agreement with pre- 
dictions made by South 35 which define the pointed cone 
half angle for which the sonift line first moves to the 
corner as a function of 7 and M m in a perfect gas at 
a = 0. (For pointed cone half angle of 70° the sonic 
line sits over the shoulder for 7 > 1.09 and M m > 10. 
For pointed cone half angle of 68° the sonic line sits 
over the shoulder for 7 > 1.12 and M m > 10.) 

The major differences between the HALIS and LAURA 
predictions are that no secondary subsonic bubbles are 
observed with HALIS and the pitching moment coeffi- 
cient is negative in all cases. (It may be that a constant 
7 solution cannot sustain the results observed in the re- 
acting gas case.) 

It should be emphasized that, although Mach num- 
ber has been used as a traditionally convenient, defin- 
ing parameter in these cases, the gas chemistry effect 
on aerodynamics is more appropriately defined in terms 
of freestream velocity. For hypersonic flows, the post- 
shock conditions over blunt bodies are a function of 
velocity and density in the freestream, and are nearly 
independent of temperature for the range of conditions 
typically encountered in atmospheric flight. Conse- 
quently, the use of Mach number introduces an inap- 
propriate freestream temperature dependence. 


(b) Windside - upper portion of sonic line is embedded in 
the captured shock. 

Figure 7: Sonic line location in the symmetry plane 
above the shoulder and cone fustrum for three Mach 
numbers spanning the static instability at a = 2° in 
the predominantly equilibrium gas chemistry regime. 
Flow enters from the top of the figure. 


Baseflow Effects on Aerodynamics 

The baseflow computation is not time accu- 
rate, but relaxation with small, constant time step 
At = .OOlj’ma.p/Ico yields time like variation. Base- 
flow computations at M m = 22.3 show no signifi- 

cant effect on aerodynamics where base pressure levels 
are approximately 0.6% of forebody levels. (Pressure 
across the base varied from 2.5 to 4.0 times freestream 
pressure.) These levels translate into a decrement of 
0.0064 in Ca at a = 5°. However, significant influence 
is evident at M m = 1.9 in Fig. 3. Unsteady recircu- 
lation patterns cause a variation in aerodynamic coeffi- 
cients bounded by the bracketed label in thesig figures. 
In this case, base pressure levels are approximately 10% 
of forebody pressure levels. 
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Aerodynamics in Nonequilibrium Flow Regime 


The effective ratio of specific heats 7 increases with 
decreasing velocity, causing a transition in the sonic 
line location across the leeside fustrum at small angles 
of attack in an equilibrium flow regime for the Mars 
Pathfinder probe. This transition may also be expected 
to occur earlier in the trajectory as the post shock gas 
chemistry changes from frozen to equilibrium and 7 
decreases. In this range, the equilibrium post shock 
density ratio is inappropriate for defining an effective 
7; consequently, there are no perfect gas computations 
offered. Lift, drag, and pitching moment coefficients, 
including laminar shear forces, are plotted as a function 
of Mach number in Fig. 8 over the hypersonic range be- 
ginning 18 seconds prior to peak heating (M m = 39) 
and ending at M m = 9.4. These results show both a 
slight decrease in axial force coefficient and a region of 
positive pitching moment coefficient at small a for the 
nonequilibrium flow domain prior to the peak dynamic 
pressure point. The trends are consistent with sonic 
line transition discussed previously. The decrease in 
axial force is associated with the pressure drop at the 
shoulder feeding upstream in the large subsonic zone. 
All of the cases at a = 2° with positive pitching mo- 
ment exhibit the leeside sonic line pattern that includes 
the secondary, subsonic bubble. This bubble sits above 
the primary sonic line that stretches along the bound- 
ary layer edge and attaches to the bow shock in front 
of the nose. 


Comparisons to Experimental Data 

The Pathfinder probe forebody shape is very sim- 
ilar to the Viking Probe shape. Only the shoulder 
radii differ (r sh, Pathfinder — 0 .02544 Pathfinder a.lld 
r sh,v iking = WlDv iking)- Although there are no 
experimental data available for aerodynamics of the 
Pathfinder shape, there are ballistic range data for the 
Viking shape obtained in COo. Data sets for C'd and 
C m ,a presented in Figs. 9 and 10 are taken from Intri- 
eri, et. at. 4 Conditions for this test are VA, k, 3.4km/s 
and Reo t co ~ 0.8 x 10 b . These conditions are very 
clos6; to the M m = 16.0 and 14.0 points in Table 1 
for Pathfinder. Comparisons in Fig. 9 show that the 
predicted values for the Pathfinder C'd are within ap- 
proximately 1%) of experimental data for the Viking 
shape for q- < 10°. Comparisons in Fig. 10 show that 
in this velocity range both experiment and calculation 
indicate a positive moment coefficient at small angles 
of attack. (No correction is made here for the different 
center-of-gravity reference point for the Viking model.) 

A calculation for the Viking shape corresponding 
to the experimental conditions of Intrieri, gt. al. 4 was 
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Figure 8: Aerodynamic coefficients versus Mach num- 
ber at several angles of attack a for both hypersonic 
nonequilibrium and equilibrium flow regimes. 
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Figure 9: Comparison of computed drag coefficient for 
Mars Pathfinder probe with measured value for geomet- 
rically similar Viking probe obtained in ballistic range 
with COo. 
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Figure 10: Comparison of computed pitching moment 
coefficient for Mars Pathfinder probe with measured 
value for geometrically similar Viking probe obtained 
in ballistic range with CCU. 


made to better validate the Pathfinder results, remov- 
ing the small differences associated with geometry, ve- 
locity, Reynolds number, and center of gravity (eg). 
The results are shown as open circles in Fig. 9 for 
Cd and hi Fig. 10 for C mQ .. As noted earlier, all 
computations assume average base pressure equal to 
Pco . The effect of higher base pressures will tend to 
lower C D (A C D = -0.02 for p hase = 3pco)- The 
pitching moment derivative is in good agreement with 
ballistic range data at small angles of attack but is ap- 
proximately 50% more negative at a = 11.5°. Rea- 
sons for this discrepancy are not understood, particu- 
larly since LAURA computations of C mQ . = —0.213 
and Ca = 1.617 at a = 11.5° are in good agree- 
ment. with flight data 3 —0.21 < C' mQ . < —0.17 and 
1.61 < Ca < 1.68 at conditions approximating peak 
dynamic pressure ( 1% = 3.4 km/s and p^ = 0.00098 


F igure 1 1 : Total angle of attack as a function of velocity 
along the trajectory as determined from a six-degree- 
of-freedom trajectory analysis using aerodynamic coef- 
ficients computed by LAURA. 


kg/rn 3 ). (The specified eg location is 0.23 diameters be- 
hind the nose and 0.01326 diameters off axis for flight 
conditions.) The sensitivity of pitching moment co- 
efficient to velocity (which affects sonic line location 
through the variation of j) for Pathfinder at a < 5° 
is clearly evident Fig. 10. 


Six- Degree-of- Freedom Trajectory 


The total angle-of-attack ax (root- mean-square 
of a and f3) as a function of velocity for the Mars 
Pathfinder probe is presented in Fig. 11 for a represen- 
tative entry. This data was produced through integra- 
tion of the six degree-of-freedom equations of motion 
with the Program to Optimize Simulated Trajectories 
using the LAURA aerodynamics and a Ma.rs-GRAM 
atmosphere. 3b The probe begins its descent at 125 km 
with a relative velocity of 7.6 km/sec and a flight path 
angle of -14.8 deg. For this trajectory, a ballistic coef- 
ficient of 60 kg/rn 2 is assumed. Entry is initiated with 
zero degree total angle-of-attack and a roll rate of 2 
rpm. Because the center-of-gravity is assumed to be 
on the symmetry axis, a nominal trim angle of zero 
degrees is expected. However, as shown in Fig. 11, 
throughout the two instability regions, significant atti- 
tude dispersions occur (with mis values of 3-4 deg). A 
more detailed examination of the Mars Pathfinder six 
degree-of-freedom flight dynamics using LAURA aero- 
dynamics has been presented. 3 ' 
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x, m 


Figure 12: Surface heating distributions on symmetry 
plane for several angles of attack a at the peak dynamic 
pressure point M m = 22.3. The windside domain is 
defined by x < 0.. 

Convective Heating 

Convective heating distributions along the symme- 
try plane at M m = 22.3 (peak dynamic pressure tra- 
jectory point) are presented in Fig. 12. Simple bound- 
ary conditions include wall temperature set to 2000 Iv 
and species recombine to freestream mass fractions at 
the wall. There is no ablation or turbulence modelling. 
The maximum momentum thickness Reynolds number 
computed from the laminar profile is approximately 575 
at peak dynamic pressure, occuring just upstream of 
the shoulder. 

Note? that both windside and leeside heating lev- 
els decrease with increasing angle of attack. At 
Mco = 22.3 and a = 0° the sonic line runs along 
the boundary layer edge, with a very slender subsonic 
bubble sitting behind the captured shock above the fus- 
trurn. The predicted stagnation point heating levels at 
o- = 0 are in good agreement (approximately 8%) 14 
with the Sutton and Graves correlation 38 for sphere 
in equilibrium flow using nose radius as a scale factor. 
This good agreement indicates that stagnation point 
velocity gradients are only weakly influenced by the 
presence of the cone. At larger angles of attack, the 
presence of the cone on the windside appears to have 
stronger influence on stagnation point gradients, lower- 
ing the heating levels because the effective body radius 
increases. 

Even though heating decreases with increasing an- 
gle of attack, it should be noted that drag coefficient 
also decreases. The probe does not decelerate as ef- 
fectively. Consequently, the peak heating and dynamic 
pressure points will occur later in the trajectory, pos- 
sibly resulting in more severe aerothermal loads on the 
vehicle if it flies at non-zero angle of attack. 


Numerical Issues 

The first, two blocks of Table 4 show the effect of 
increasing surface grid resolution by 75%). The change 
in Ca is less than 0.25%). Changes in C/v and C m are 
16%) and 33%) respectively for the M m = 16 point at 
a = 2°. These values are more sensitive to grid be- 
cause they represent the difference of two large, nearly 
equal contributions from the windside and leeside of 
the probe. Note, however, that these changes are rela- 
tively small (2.9%) and 6.3%), respectively) when mea- 
sured relative to the change in value going from a = 2° 
to q- = 5°. 

The contribution of viscous forces to the aerody- 
namics at small angles of attack are very small over 
most of thgTrajectory, as would be expected for a blunt 
body. The integration of only the pressure force contri- 
bution to aerodynamics is presented in the second and 
third blocks of Table 4 for two representative points. 

Differences between the LAURA and HALIS re- 
sults are predominantly due to the gas model. This 
sensitivity may be confirmed by examining the fourth 
and fifth blocks of Table 4. If LAURA is run in an 
inviscid, perfect gas mode like HALIS, then the differ- 
ences in solution are algorithmic in nature. For the two 
representative points tested at a = 5°, differences in 
sensitive quantities like C m for the two perfect gas so- 
lutions are less than 4%) at M m = 22.3 and less than 
10%) at Mco = 16.0. When the real gas model and 
viscous terms are included in LAURA, the effects on 
aerodynamics are much larger (24%) at M m = 22.3 
and 25%) at M m = 16.0). Since viscous contributions 
have already been shown to be negligible, the sensitiv- 
ity of aerodynamics on the gas model is established. 


Conclusions 

(1) As the Mars Pathfinder probe descends 
through the Martian atmosphere the minimum value 
of the post-shock effective 7 first, decreases from frozen 
gas chemistry values (ps 1.333) to equilibrium values 
(1.094) corresponding to a. velocity of 4.86 krn/s. As 
the probe continues to decelerate through an equilib- 
rium post-shock gas chemistry regime, the value of 7 
increases again, until reaching its perfect, gas value of 
1.333 at parachute deployment. (0.42 krn/s). 

(2) At. small angles of attack (o- < 5°) the sonic 
line location shifts from the shoulder to the nose cap 
and back again on the leeside symmetry plane because 
of the change in 7 and the cone half-angle of 70°. A 
noteworthy aspect, of this transition process is the for- 
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mation of a secondary, subsonic zone over the leeside 
fustrum. 

(3) Pressure distributions on the cone fustrum ap- 
proaching the shoulder tend to be very flat when the 
sonic line sits forward over the spherical nose. Effects 
of the expansion over the shoulder can only be commu- 
nicated upstream through the subsonic portion of the 
boundary layer. In contrast, pressure distributions on 
the cone fustrum approaching the shoulder tend to be 
more rounded when the sonic line sits on the shoulder, 
exhibiting a more pronounced influence of the expan- 
sion on the upstream flow. 

(4) In general, windside pressures exceed leeside 
pressures on the cone fustrum, producing a stabilizing 
moment that pitches the probe back to zero angle of 
attack. However, the behavior of the pressure distri- 
bution in the vicinity of the shoulder significantly in- 
fluences the pitching moment coefficient because of the 
relatively larger moment arms and surface area at the 
edge of the probe as compared to the inboard nose and 
fustrum regions. For the Mars Pathfinder probe at 2° 
angle of attack, the flat, leeside pressures approaching 
the shoulder (when the sonic line sits over the nose) can 
exceed the rounded windside pressures approaching the 
shoulder (when the sonic line sits over the shoulder). 
The net effect of this crossover distribution near the 
shoulder tends to pitch the probe to higher angles of 
attack. The overall balance (crossover-point location) 
is sensitive to both freestream conditions and the gas 
chemistry. The combination of small angle of attack, 
large cone angle, windside sonic line over the shoulder, 
and leeside sonic line over the nose with a secondary, 
subsonic zone above the leeside fustrum yield condi- 
tions which favor (but do not necessarily guarantee) a 
destabilizing pitching moment. 

(5) Conditions for a positive, destabilizing moment 
coefficient derivative C mt0l occur twice in the Mars 
Pathfinder mission as determined by the most accurate 
flow simulation used in this study. The first occurence 
(7.5 > VC, > 6.5 km/s, 51 > h > 37 km, vicinity 
of peak heating for this trajectory) results from the 
transition in the sonic line location as a function of 
gas chemistry changing from nonequilibrium to equi- 
librium. The second occurence (4.0 > Coo > 3.1 km/s, 
25 > h > 22 km) results from the transition in the sonic 
line location as a function of decreasing flow enthalpy 
in an equilibrium gas chemistry regime. 

(6) Sonic line movement affects the heating dis- 
tributions by altering the effective radius of curvature 
of the body. Peak heating tends to decrease as angle 
of attack increases for fixed freestream conditions on 
the Mars Pathfinder probe because the body appears 
blunter to the oncoming flow. However, the drag co- 
efficient decreases with increasing angle of attack, so 
that the ballistic coefficient increases and peak heating 
point along the trajectory may be more severe. 


(7) A sample, six-degree-of-freedom trajectory 
based on the computed aerodynamics and a spin rate of 
2 revolutions per minute predicts a precessing motion 
during the two, statically unstable periods in which the 
root-mean-square value of angle of attack is between 3 
and 4 degrees. 
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Table 1 - Trajectory Points 


Goo rn/s 

Poo , kg/m 3 

Too, K 

Twall, K 

h, km 

Moo 

■.JJMMliMI 


7 

comment 

7623. 

9.41 10 -6 

139.5 

1600. 

69.7 

39.3 

0.017 

- 

- 

- 

7407. 

9.48 10 -5 

151.7 

1600. 

50.6 

36.7 

0.158 

- 

- 

- 

6592. 

3.24 10 -4 

162.0 

2000. 

40.7 

31.6 

0.460 

- 

- 

max. htg. 

5596. 

7.91 10 -4 

152.2 

2000. 

31.4 

27.7 

0.992 

- 

- 

- 

4862. 

1.16 10 -3 

176.8 

2000. 

28.5 

22.3 

1.149 

21.4 

1.0938 

max. Q 

4521. 

1.28 10 -3 

179.0 

2000. 

26.7 

20.9 

1.168 

- 

- 

- 

3996. 

1.55 10 -3 

181.6 

2000. 

24.8 

18.3 

1.239 

19.7 

1.1010 

- 

3515. 

1.85 10 -3 

183.8 

2000. 

23.0 

16.0 

1.290 

17.9 

1.1105 

- 

3085. 

2.14 10 -3 

186.0 

2000. 

21.5 

14.0 

1.300 

- 

- 

- 

2707. 

2.44 10 -3 

188.4 

2000. 

20.2 

12.2 

1.290 

13.9 

1.1405 

- 

2096. 

3.03 10 -3 

191.5 

2000. 

18.0 

9.4 

1.227 

10.8 

1.1800 

- 

- 

- 

- 

- 

- 

6.6 

- 

8.5 

1.2125 

- 

- 

- 

- 

- 

- 

2.0 

- 

2.8 

1.3300 

- 

420. 

7.49 10 -3 

200.0 

300. 

8.2 

1.9 

0.591 

- 


parachute 



Table 2 - Perfect Gas Aerodynamic Coefficients 


a, deg 

Moo 

Cn 

C A 

C m 

c D 

c L 

0. 

22.3 

0.00000 

1.7160 

0.00000 

1.7160 

0.00000 

0. 

18.3 

0.00000 

1.7150 

0.00000 

1.7150 

0.00000 

0. 

16.0 

0.00000 

1.7122 

0.00000 

1.7122 

0.00000 

0. 

12.2 

0.00000 

1.6803 

0.00000 

1.6803 

0.00000 

0. 

9.4 

0.00000 

1.6358 

0.00000 

1.6358 

0.00000 

0. 

6.6 

0.00000 

1.5970 

0.00000 

1.5970 

0.00000 

0. 

2.0 

0.00000 

1.3426 

0.00000 

1.3426 

0.00000 

2. 

22.3 

0.00678 

1.7100 

-0.00431 

1.7091 

-0.05290 

2. 

18.3 

0.00563 

1.7090 

-0.00316 

1.7081 

-0.05401 

2. 

16.0 

0.00427 

1.7040 

-0.00179 

1.7031 

-0.05520 

2. 

12.2 

0.00297 

1.6767 

-0.00093 

1.6757 

-0.05555 

2. 

9.4 

0.00421 

1.6347 

-0.00255 

1.6338 

-0.05284 

2. 

6.6 

0.00482 

1.5960 

-0.00322 

1.5951 

-0.05088 

2. 

2.0 

0.00458 

1.3421 

-0.00338 

1.3414 

-0.04226 

5. 

22.3 

0.01573 

1.6836 

-0.00929 

1.6785 

-0.13106 

5. 

18.3 

0.01433 

1.6809 

-0.00804 

1.6757 

-0.13222 

5. 

16.0 

0.01280 

1.6764 

-0.00674 

1.6711 

-0.13335 

5. 

12.2 

0.00922 

1.6618 

-0.00404 

1.6562 

-0.13565 

5. 

9.4 

0.01013 

1.6278 

-0.00599 

1.6224 

-0.13178 

5. 

6.6 

0.01179 

1.5890 

-0.00792 

1.5839 

-0.12674 

5. 

2.0 

0.01147 

1.3370 

-0.00845 

1.3329 

-0.10510 

8. 

22.3 

0.02587 

1.6260 

-0.01572 

1.6137 

-0.20067 

8. 

18.3 

0.02467 

1.6428 

-0.01476 

1.6302 

-0.20420 

8. 

16.0 

0.02323 

1.6397 

-0.01356 

1.6269 

-0.20519 

8. 

12.2 

0.01975 

1.6265 

-0.01117 

1.6134 

-0.20680 

8. 

9.4 

0.01748 

1.6059 

-0.01031 

1.5927 

-0.20618 

8. 

6.6 

0.01852 

1.5740 

-0.01226 

1.5612 

-0.20071 

8. 

2.0 

0.01904 

1.3289 

-0.01383 

1.3186 

-0.16609 

11. 

22.3 

0.03738 

1.6010 

-0.02341 

1.5787 

-0.26879 

11. 

18.3 

0.03647 

1.5970 

-0.02268 

1.5746 

-0.26892 

11. 

16.0 

0.03520 

1.5935 

-0.02192 

1.5709 

-0.26950 

11. 

12.2 

0.03198 

1.5811 

-0.01984 

1.5581 

-0.27029 

11. 

9.4 

0.02855 

1.5657 

-0.01787 

1.5423 

-0.27072 

11. 

6.6 

0.02652 

1.5460 

-0.01733 

1.5226 

-0.26895 

11. 

2.0 

0.02563 

1.3144 

-0.01873 

1.2951 

-0.22564 
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Table 3 - Real Gas Aerodynamic Coefficients 


a, deg 

Moo 

Cn 

C A 

C m 

c D 

C L 

0 . 

31.6 

0.00000 

1.6984 

0.00000 

1.6984 

0.00000 

0 . 

22.3 

0.00000 

1.7290 

0.00000 

1.7290 

0.00000 

0 . 

9.4 

0.00000 

1.6393 

0.00000 

1.6393 

0.00000 

0 . 

1.9 

0.00000 

1.3079 

0.00000 

1.3079 

0.00000 

2. 

39.3 

0.00749 

1.6597 

-0.00122 

1.6589 

-0.05044 

2. 

36.7 

0.00358 

1.6796 

0.00005 

1.6787 

-0.05504 

2. 

31.6 

0.00238 

1.6943 

0.00110 

1.6934 

-0.05675 

2. 

27.7 

0.00623 

1.7138 

-0.00344 

1.7130 

-0.05358 

2. 

22.3 

0.00621 

1.7208 

-0.00369 

1.7199 

-0.05385 

2. 

20.9 

0.00557 

1.7168 

-0.00317 

1.7159 

-0.05435 

2. 

18.3 

0.00344 

1.7141 

-0.00110 

1.7131 

-0.05639 

2. 

16.0 

0.00173 

1.6999 

0.00066 

1.6989 

-0.05760 

2. 

14.0 

0.00262 

1.6740 

-0.00080 

1.6730 

-0.05580 

2. 

12.2 

0.00322 

1.6547 

-0.00160 

1.6538 

-0.05453 

2. 

9.4 

0.00364 

1.6382 

-0.00218 

1.6373 

-0.05353 

2. 

1.9 

0.00480 

1.3069 

-0.00342 

1.3062 

-0.04081 

5. 

39.3 

0.01852 

1.6508 

-0.00273 

1.6461 

-0.12543 

5. 

36.7 

0.01121 

1.6595 

-0.00290 

1.6542 

-0.13347 

5. 

31.6 

0.01047 

1.6616 

-0.00366 

1.6562 

-0.13439 

5. 

27.7 

0.01451 

1.6821 

-0.00723 

1.6769 

-0.13215 

5. 

22.3 

0.01425 

1.6877 

-0.00735 

1.6825 

-0.13289 

5. 

18.3 

0.01104 

1.6804 

-0.00576 

1.6749 

-0.13546 

5. 

16.0 

0.00996 

1.6679 

-0.00460 

1.6624 

-0.13545 

5. 

14.0 

0.00859 

1.6563 

-0.00414 

1.6507 

-0.13580 

5. 

12.2 

0.00890 

1.6418 

-0.00466 

1.6363 

-0.13422 

5. 

9.4 

0.00971 

1.6278 

-0.00569 

1.6224 

-0.13220 

5. 

1.9 

0.01190 

1.3027 

-0.00846 

1.2987 

-0.10169 

11. 

39.3 

0.04517 

1.5896 

-0.01247 

1.5690 

-0.25897 

11. 

36.7 

0.03592 

1.5767 

-0.01797 

1.5546 

-0.26558 

11. 

22.3 

0.03571 

1.5954 

-0.02202 

1.5729 

-0.26937 

11. 

1.9 

0.02607 

1.2813 

-0.01850 

1.2627 

-0.21890 


Table 4 - Aerodynamic Coefficients 


Block 

Code 

type 

grid 

a, deg 

Moo 

Cn 

C A 

C m 


LAURA 

real vis. 

61x31x65 

0 . 

22.3 

0.00000 

1.7290 

0.00000 

1 

LAURA 

real vis. 

81x41x65 

0 . 

22.3 

0.00000 

1.7290 

0.00000 


LAURA 

real vis. 

61x31x65 

2. 

16.0 

0.00173 

1.6999 

0.00066 


LAURA 

real vis. 

81x41x65 

2. 

16.0 

0.00149 

1.7036 

0.00100 

2 

LAURA 

pressure 

81x41x65 

2. 

16.0 

0.00150 

1.7016 

0.00088 


LAURA 

real vis. 

61x31x65 

0 . 

31.6 

0.00000 


.6984 

0.00000 

3 

LAURA 

pressure 

61x31x65 

0 . 

31.6 

0.00000 


.6944 

0.00000 


LAURA 

real vis. 

61x31x65 

5. 

22.3 

0.01425 


.6877 

-0.00735 


LAURA 

perf. inv. 

61x31x33 

5. 

22.3 

0.01523 


.6894 

-0.00963 

4 

HALIS 

perf. inv. 

46x37x21 

5. 

22.3 

0.01573 


.6836 

-0.00929 


LAURA 

real vis. 

61x31x65 

5. 

16.0 

0.00996 


.6679 

-0.00460 


LAURA 

perf. inv. 

61x31x33 

5. 

16.0 

0.01155 


.6788 

-0.00616 

5 

HALIS 

perf. inv. 

46x37x21 

5. 

16.0 

0.01280 


.6764 

-0.00674 


14 









































































